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Abstract. Due to the continuous developments in polarimetric instrumentation, 
which will become even more dramatic in the near future with the availability of new 
generation solar telescopes, we are now severely confronted with a variety of new de- 
tailed observations of high diagnostic potential, whose interpretation requires a firmly 
established theoretical framework. In this contribution, I review the fundamental phys- 
ical processes that underlie the generation and transfer of polarized radiation in stellar 
atmospheres, and I discuss the present status of the theoretical schemes now available, 
pointing out their main successes and limitations. I also present some ideas about the 
theoretical improvements that I consider necessary to achieve a correct interpretation of 
the complex phenomenology shown by polarimetric observations, focusing particularly 
on the second solar spectrum, which can be considered as one of the most important 
test benches of the theory. 



1. Introduction 

Polarized radiation is expected whenever the interaction process between matter and 
radiatio n is characterized by any kind of sym metry-breaking (see, for example, the re- 
view by Casin i & Landi Deal ' In nocent j|2007h . The symmetry-breaking can be intrinsic 
to the interaction itself, or due to external causes (e.g., presence of a magnetic field). 
In order to get more insights on this fundamental statement, let us consider a quantum 
transition between an upper level J u , and a lower level J[, J being the quantum number 
associated with the total angular momentum operator. According to quantum mechan- 
ics, the atomic system in the upper (or lower) level can occupy any of the (27 + 1) 
possible M magnetic substates, which are degenerate if no magnetic field is present. 
The radiation emitted by the various transitions between the M u substates of the upper 
level and the M[ substates of the lower level (generally referred to as Zeeman com- 
ponents), is characterized by well defined polarization properties which depend on the 
angle between the emission direction and the quantization axis, and on the value of 
AM - M u - M[. In the electric dipole approximation the only allowed transitions are 
those with AM = (jt components), and those with AM = ±1 (cr\, and <x r components, 
respectively). The polarization properties of the n and cr components emitted along 
various directions are schematically shown in Fig. 1. 

In a de-excitation process, each Zeeman component contributes to the emitted 
radiation with a particular weight which takes into account the relative population of 
the upper M-sublevel of the component, the relative strength of the component, and the 
emission direction. If the atomic system is naturally (isotropically) excited, then the 
upper magnetic sublevels are equally populated, and it can be shown that the weights 
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Figure 1 . Polarization properties of the radiation emitted by the different Zeeman 
components. The z axis is the quantization axis. 



of the various Zeeman components assume values such that, if no magnetic field is 
present, the total radiation emitted in any direction is completely unpolarized. 

There are two basic scenarios (not mutually exclusive) where the different polar- 
ization properties of the Zeeman components manifest themselves in emission of po- 
larized radiation. The first scenario corresponds to the possibility that the M substates 
may be separated in energy, so that the contributions of the various components fall at 
slightly different wavelengths. In this case, even if the the atomic system is naturally 
excited, the emitted radiation has polarization properties varying with wavelength. The 
most common circumstance for this condition to occur is when the atomic system inter- 
acts with an external magnetic field {Zeeman effect). The second scenario corresponds 
to the possibility that the upper M substates are unevenly populated, so that the various 
components no longer contribute to the emitted radiation with those particular weights 
that make the total polarization vanish. This condition generally occurs whenever the 
atomic system is excited through a physical process which, for any reason, is not spa- 
tially isotropic. An atom whose magnetic sublevels are not evenly populated and/or are 
characterized by well-defined phase relations is said to be polarized. The radiation emit- 
ted by a polarized atomic system is in general polarized, as it can be shown through the 
following clarifying example. Let us consider a two level atom, the lower level having 
J i - and the upper level having J u = 1 , and let us suppose to excite it with a collimated 
unpolarized radiation beam propagating along the quantization axis. In the presence of 
such an incident radiation field, it can be shown that only the M u = ± 1 magnetic sub- 
levels will be populated, so that only the cr components will contribute to the emitted 
radiation. From Fig. 1 it is then clear that if the radiation emitted perpendicularly to the 
incident beam is considered, this will be totally linearly polarized perpendicularly to 
the scattering plane. It should be observed that this is basically the physical mechanism 
responsible for the linearly polarized spectrum of the solar radiation coming from quiet 
regions close to the limb (the so-called Second Solar Spectrum). Of course, in the solar 
atmosphere the radiation field impinging on the atoms is not perfectly collimated along 
the local vertical, nevertheless it has a certain degree of anisotropy that provides for the 
radiation scattered at 90° a linear polarization degree of the order of 1%. Polarization 
produced through resonance scattering processes like the one here described is usually 
referred to as scattering polarization. 

The quantum theory of the Zeeman effect was developed in the first half of the last 
century, and we can safely say that today it is absolutely well established. Indeed, the 
difficulties related to the interpretation of spectropolarimetric signals produced through 
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this mechanism are essentially due to the complexity of the line formation problem in 
the solar atmosphere, the physical origin of the polarization being, on the other hand, 
absolutely clear. On the contrary, as soon as the first scattering p olarization signals, 
detected at the solar limb since the early 1940s (see iRedmanll 194-lf) . could be observed 
with a sufficient accuracy, it was immediately clear that the theoretical interpretation 
of the peculiar profiles shown by several of them would have required a substantial 
improvement of our understanding of the physics of resonance scattering. The incom- 
plete comprehension of this fundamental physical process has become particularly ev- 
ident during t he last deca des when, thanks to the development of polarimeters like 
ZIMPOL (see lPovelll995 ). able to reach sensitivities of the order of 1 part over 10 5 , 
many new app arently "eni gmatic" scattering polarization signals have been detected 
(see lStenflo & Kellerll 19971) . Given the high number of physical "ingredients" involved 
in the generation of scattering polarization signals, and given their consequent enor- 
mous diagnostic content, since the 1970s their interpretation has been representing one 
of the most important and exciting theoretical challenges in the field of solar polarime- 
try. 

Though our understanding of the physics of resonance scattering has substantially 
improved during the last 30 years, and a robust theoretica l scheme, based on the princi- 
ples o f quantum electrodynamics has been developed (see lLandi DeglTnnocenti & Landolfi 
2004, hereafte r LL04), and success fully applied for the interpretation of a wide range 
of signals (see Trujillo Buenol l2009l for a recent review), our comprehension of the rich 
phenomenology shown by scattering polarization signals remains rather fragmentary. 

In the present contribution I will review some of the most important theoretical 
approaches to the physics of scattering polarization developed so far. The aim of the 
paper is to point out strengths and limitations of the various approaches by reviewing 
the physical processes and phenomena they are able to describe, and by discussing 
the most important hypotheses and approximations they are based on (for a detailed 
derivation of each approach the reader will be referred to the original papers or to more 
exhaustive monographs). Examples of the main successes of each theory, as well as a 
discussion of unsolved problems, will also be presented. Of course, the review is by no 
means exhaustive, and I apologize in advance to all whose work is not mentioned here. 



2. The second solar spectrum 

The linearly polarized solar limb spectrum, today known as the second solar spectrum, 
is probably the clearest manifestation of scattering polarization in astrophysics, and it is 
no coincidence that it played, and is still playing a fundamental role in the development 
of the theory of polarization. Its importance can be clearly highlighted through a few 
remarkable examples of signals whose interpretation required fundamental advances in 
our understanding of the physics of resonance scattering. 

The interpretation of the peculiar polarization pattern observed across the H and 
K lines of Ca ii (and, similarly, across the Di and D2 lines of Na 1), which shows a sign 
reversal between the two lines, required the comprehension of the fundamental role 
played by quan tum interferenc es (or coherences) between different quantum levels, as 
first showed by IStenflol (Il980h . In 1997, Stenflo explained the triplet-peak structure of 
the second solar spectrum of the Ba 11 D2 line in terms of hyperfine structure (HFS), 
pointing out how this physical aspect, generally negligible in the usual intensity spec- 
tram, can be extremely important when polarization phenomena are taken into account 
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(see IStenflolll997l) . In the same year, iTrujillo Bueno & L andi Pegl'In nocenti! (1 19971) 



pointed out the importance of lower level polarization, previously systematically ne- 
glected due to the qualitative idea that atomic polarization in long-lived levels should be 
always completely destroyed by depolarizing collisions and magnetic fields. By taking 
into account this new ingredient it has been possible to i nterpret the polarization signals 
obser ved in the Ca n infrared triplet line at 8662 A (see iManso Sainz & Trujillo Buenol 



2003), as well as in the Mg i b-lines (see lTrujillo Buenolll999l 12001). 



There are now clear evidences that the interpretation of a large class of signals (like 
the one observed in the Ca i line at 4227 A) requires the effects of partial redistribution 
in frequency to be taken into account. The development of theoretical frameworks able 
to account for these effects is nowadays one of the hottest topics in this research field. 
Laboratory experiments are now being performed in order to clarify various aspects of 
this spectrum, like the pres ence of line-integrated line ar polarization in the Di line of 
sodium and potassium (see Thalman n et alj|2006ll2009l) . and in order to verify the main 
results derived from the theoretical approaches developed so far. 

The importance of the second solar spectrum is not limited to s ingle signals, but 



d2i 

120051) . iBelluzzi & Landi Pegl'Innocenti! d2009l) identified the strongest signals of this 
spectrum, and divided them into five different classes according to the shape of their 
profile. Through such a systematic analysis, the authors discovered a series of spectro- 
scopic properties common to the spectral lines responsible for such signals, properties 
that have been summarized into the following three empirical laws: 

1. First law: the transitions producing the strongest polarization signals of the sec- 
ond solar spectrum are either resonance transitions, or subordinate transitions 
whose lower level is the upper level of a resonance transition producing a strong 
polarization signal. 

2. Second law: all the strong polarization signals that are produced by spectral 
lines having a small equivalent width (i.e. lines having Wa/A < 20 F) are of type 

3. Third law: the spectral lines producing the strongest polarization signals ( Q/I > 
0.17%) are due to transitions having either AJ — J u — J( = +1 or AJ = 0. 

The theoretical interpretation of these laws is already under investigation, and it will 
probably represent a further important test bench for the theoretical approaches that 
will be developed in the future. 



also concern s some of its general properties, as recently pointed out by Belluzzi & Landi Pegl'Innocenti 
(120091 . | 2010l>. Using the atlas "The Second So lar Spectrum" dGandorferll2000L 12002 . 



3. A unifying theoretical approach in the formalism of quantum electrodynamics 

The unifyin g, self-consistent theo r etical approach to the physics of polarization, for- 
mulated by Landi Pegl'Innocenti! dl983l ) within the framework of the density-matrix 
formalism, starting from the principles of quantum electrodynamics, represents one of 
the most important theoretical accomplishments obtained so far in this research field. 



'All the signals whose profile shows a single sharp peak have been classified as "S" signals. 
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A clear and detailed derivation of this theoretical scheme can be found in LL04, and 
will not be repeated here. In this section, in order to point out strengths and weaknesses 
of this approach, I will present only the fundamental equations the theory is based on, 
discussing their physical meaning, and focusing the attenti on on the most important hy- 



jotheses and approximation s that are introduced (see also lCasini & Landi Degl'Innocentil 
120071 : iTrujiilo Buenolll990h . I will then review some of the most important successes 
of this theory, and I will discuss its main limitations and their consequences for the 
interpretation of particular classes of signals. 

Other theoretical approaches, developed with different assumptions and suitable 
for investigation of those signals which cannot be interpreted under the approximations 
of the scheme described here, will be presented in the next section. 



3.1. The density operator 

The first problem to be addressed in developing a quantum theory for scattering po- 
larization is to find a suitable way to describe atomic polarization (i.e. the presence 
of population unbalances and/or of quantum interfere nces between pairs of m a gnetic 
sublevel s). A very good solution, as proposed by iBommier & Sahal-Brechotj dl978h 
and by IBommier] (Il980l) is through the density operator, an extremely powerful tool 
for describing any physical system that is in a statistical mixture of states. The density 
operator is defined by 

P = ^Pa\^a)(^a\ , (1) 



where p a is the probability for the system to be in the (pure) dynamical state described 
by the vector \\j/ a ), and where the sum is extended to all the pure states in which the 
system can be found. The action of the density operator is completely specified once 
its matrix elements, evaluated on a given basis of the Hilbert space associated to the 
quantum system, are known. Such density-matrix elements contain all the accessible 
information about the system; through them it is possible to describe the dynamical 
state of the system and its temporal evolution in a very compact way. For an atomic 
system, the most natural basis for defining the matrix elements of the density operator is 
the basis of the eigenvectors of the total angular momentum. On this basis the elements 
of the density matrix are given by 

(aJM I p I a'J'M') = p (aJM, a'J'M') , (2) 

where / and M are the angular momentum quantum numbers previously introduced, 
while a represents a set of inner quantum numbers. As it can be shown (see, e.g., LL04), 
the diagonal elements represent the populations of the magnetic sublevels, while the 
off-diagonal elements describe the quantum interferences between different magnetic 
sublevels. The expectation value of a given dynamical variable O is given by (see, e.g., 
LL04) 

(0)= Tr{p 6), (3) 



with O the quantum operator associated to the dynamical variable O. 
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3.2. The statistical equilibrium and the radiative transfer equations 

The time evolution of the density operator associated to a physical system described by 
the Hamiltonian H is given by the Liouville equation 

d 27ri r 

If the physical system is composed by one atom and by the radiation field, then the 
Hamiltonian H is the given by the sum H = Ha + Hr + V, with Ha the atomic Hamil- 
tonian (which eventually includes the interaction with an external magnetic field), Hr 
the Hamiltonian of the radiation field (generally described within the formalism of the 
second quan tization), and V the interactio n Hamiltonian (for its explicit expression see, 
for example, ICohen-Tannoudji et all 19771) . 



When the total Hamiltonian can be expressed as the sum of an unperturbed Hamil- 
tonian and an interaction Hamiltonian, it is generally useful to work in the so-called 
interaction picture. In this representation, in fact, only the interaction Hamiltonian ap- 
pears in the commutator in the right-hand side of the Liouville equation, which assumes 
the form 

d 27ri r 

- Pl (t) = -—[V hPl (t)] , (5) 

where the label / means that the corresponding quantity is defined in the interaction 
picture. Equation © has the formal solution 

Pi(0=Pi(0)-^ rVi(0,Pi(0]df'. (6) 
h Jo 

This equation can be used again to express the quantity p\{t') that appears in the inte- 
gral. By repeating this procedure, a perturbative development of the formal solution of 
Eq. (f5]) is obtained. If the development is stopped at the fc-th term, by substituting the 
quantity pi(f') in the integral with pi(0), the matter-radiation interaction is said to be 
described at the k-th perturbative order. 

The equation describing the time evolution of the expectation value of a given 
dynamical variable O can be obtained calculating the time derivative of Eq. (f3j, and 
using Eq. Q. After some algebra one obtains 



- 0(t) = Tr{(- Oi(0)pi(0} - -^Tr{[d!(f), Vi(0]pi(0)} 

-^Tr{ J [[Oiit), Vk(t)], Vi(?')]pi(Odf'} . (7) 



This is an exact equation, and represents the starting point to derive the equations de- 
scribing the evolution of the atomic system and of the radiation field resulting from their 
mutual interaction. In particular, if in place of 0(t) and Oi(t) one puts the density ma- 
trix elements associated to the atomic system, and the corresponding quantum operator, 
then the statistical equilibrium equations for the density matrix are obtained. If in place 
of 0(t) and 0\(t) one puts the Stokes parameters describing the radiation field, and the 
corresponding quantum operators, then the radiative transfer equations are obtained. 

The approximation that is now introduced is to substitute the quantity pi(?') that 
appears in the integral in the right-hand side of Eq. (|7) with pi(0), so that it can be 
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extracted from the integral. Since this is equivalent to performing the same substitu- 
tion in the integral on the right-hand side of Eq. ©, this approximation is equivalent 
to consider a i-st order perturbative development of the matter-radiation interaction, 
which means that only i-st order processes (processes involving only one photon) can 
be described. 

Within this approximation, scattering, which is intrinsically a 2-nd order process, 
has to be treated as a succession of independent i-st order absorption and emission 
processes (i.e. in the limit of complete redistribution in frequency, CRD). It should be 
observed that this is correct either in the presence of collisions that perturb the atom 
during the interaction with the radiation field to the point of completely relaxing the 
coherence of the scattering process, or when the radiation field illuminating the atom is 
spectrally flat (i.e. independent of frequency). Since polarization phenomena are gen- 
erally negligible when the first condition occurs (because of the effect of depolarizing 
collisions), the latter (usually referred to as the flat-spectrum approximation) is gener- 
ally assumed when the discussed theoretical approach is applied. More precisely, this 
approximation requires the radiation field to be flat across spectral intervals larger than 
the Bohr frequency connecting pairs of magnetic sublevels between which quantum 
coherences are considered, and larger than the inverse lifetime of the various magnetic 
sublevels. 

The statistical equilibrium equations for th e density matrix were first derive d (un- 
der the previous ly discussed approximation) by iBommier & Sahal-Brechotl (Il978h and 
lBommierl(ll980l) . The equations describe the transfer and relaxation of populations and 
coherences due to absorption and emission processes, and due to the presence of a mag- 
netic field (Hanle effect). The same e quations were then rederived , under the same ap- 
proximation, in the seminal paper by lLandi DeglTnnocentil (119831) . In this work Landi 
DeglTnnocenti derived within the same formalism both the statistical equilibrium equa- 
tions for the density matrix, and the radiative transfer equations for polarized radiation 
(described in terms of the four Stokes parameters), thus "closing" the Non-LTE prob- 
lem in the presence of polarization phenomena (Non-LTE problem of the 2 nd kind). 
The equations include the effects of a magnetic field (Zeeman and Hanle effects), and 
of collisions (elastic and inelastic). A detailed derivation of these equations, as well as a 
complete discussion of the important aspect (here neglected) of velocity/density-matrix 
correlations, can be found in LL04. 

3.3. Strengths and successes of the theory 

The main strength of the theoretical scheme described above is clearly its self-consis- 
tency, since all the equations are derived within the same elegant formalism from the 
principles of quantum electrodynamics. Going to more specific aspects, particularly 
important strength points of the theory are the following: 

• possibility to take into account lower level polarization; 

• possibility to take into account coherences within the same /-level (multi-level 
atoms) and between different /-levels (multi-term atoms); 

• possibility to treat atoms with hyperfine structure; 

• possibility to describe the interaction with a magnetic field in every regime, from 
the Zeeman effect regime to the complete Paschen-Back effect regime; 
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• possibility to take into account both the Zeeman and the Hanle effects. 

The theory has been extensively and very successfully used for the determination of 
magnetic fields in filaments and prominences, where the flat-spectrum approximation is 
generally well justified (see , e.g..|Bommier et al.ll98 ll: lAthay et al.lll983l : lBommier et al 

)2j). 



19941 : iTrujillo Bueno et all 120021) . Thanks to the possibility of taking into account 



quantum coherences between different /-levels, through this theoretical scheme it is 
possible to reproduce the peculiar polarization patterns (see Sect. 2) observed in the 
second solar spectrum across the H and K lines of Ca n, and across the Di and D2 
line s of Na 1 (see LL04). It wa s wi thin the framework of this theoretical approach 
that lTruiillo Buenol (I1999L l200lh and iManso Sainz & Truiillo Buenol (120031) . as men- 



tioned in Sect. 2, first explained, in terms of lower level polarization, the "enigmatic" 
signals produced in the same spectrum by the Mg 1 b-lines and by the Ca 11 infrared 
triplet, respectively. The possibility of treating multi-level atoms has made this ap- 
proach particularly attractive for investigating complex atoms. Very interesting results 
hav e been obtained, for example, for the interpretation of t he second solar spectrum of 
Ti 1 dManso Sainz & Landi Degl'Innocentill2002h and Ce n dManso Sainz et alj |2006). 



An interesting result, which could be obtained thanks to the possibility for the the- 
ory to model atoms with HFS, and to take into account both the Hanle and Zeeman 
effects in every magnetic regime, is the discovery of a differential magnetic sensitivity 
of the triple-peak signal sh own in the second solar spectrum by the Ba 11 D2 line (see 
Belluzzi et al.ll2006l 120071) . Such differential magnetic sensitivity, which has great di- 



agnostic p otentialities, has already been confirmed by v arious observational campaigns 
(see, e.g. jRamelh et al.ll2009l : Edpez Ariste et al.ll2009h . It should be observed that the 



possibility of considering coherences between different /-levels, or between different 
HFS F-levels, in the presence of a magnetic field in the incomplete Paschen-Back effect 
regime allows the theory to account for q uantum mechanisms like the level-crossing and 
the anti-level-crossingP] effects (see, e.g. jBelluzzi et al]|2007l). 



Within the framework of this theory, IBeUuzzi eT al. (2009 ) investigated th e phys 



ical origin of the weak polarization signal observed by Stenflo et alj ( 2006T) in the 
second solar spectrum of the faint lithium doublet at 6708 A, and concluded, in dis- 
agreement with this observation showing a triplet-peak structure, that only a two-peak 
profile can be theoretically expected, and that only a one-peak profile should be mea- 
surable unless the temperature of the atmospheric layer the radiation comes from is 
sufficiently low. On the other hand, improved observations (like those presented in 
Stenflo 2011, and other s which are st i ll unp ublished) seem to agree with the one- 
peak p rofile predicted by IBeUuzzi et al. (2009). The theoretical result of IBeUuzzi et alj 
(2009) may thus be considered as an example of a theoretical prediction that anticipated 
an observational finding, an important indication of the adequacy of the theory. The 
soundness of the theory can be also appreciated in that it has allowed to predict new 
physical mechanisms like the so-called "alignment-to-orientation conversion mecha- 
nism" (see LL04), or the "atom ic -polarization-induced Faraday pulsation effect" (see 
lLandi DeglTnnocenti et"aT]|201ll) . 

It should be finally pointed out that the results obtained by applying this theory to 
different multiplets are in good qualitative agreement with the third law formulated by 



2 This expression has been proposed by Bommie 
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this theoretical scheme can be found in 



Belluzzi & Landi DeglTnnocentil d2009 ). Other examp l es of successful applications of 



Truii llo Buenol (l2009h . 



3.4. Limitations and open problems 

The main limitation of this theoretical approach is the impossibility to take into account 
effects of partial redistribution (PRD) in frequency. While in the intensity spectrum 
such effects are important (mainly in the wings) only for a limited number of very strong 
resonance lines (e.g., Ca n H&K, Mg n h&k, Ha, Lya), in the seco nd solar spectrum, a s 
pointed out by several works, starting with the pioneering one by lOmont et al.1 dl972l) . 
they seem to be the key ingredient to i nterpr et most of the linear polarization profiles 
that Bel luzzi & Landi Degl'Innocentil ((2009) have classified as "M" signal^j (typical 
examples are the signals produced by the Na i D2 line or by the Ca 1 line at 4227 A). 

The open problems, however, are not limited to the interpretation of the "M" sig- 
nals in terms of PRD effects. An interesting example of an unsolved problem concern- 
ing the second solar spectrum is the physical origin of the triplet-peak profile shown 
by the Sc 11 line at 4247 A. This signal is very similar to the one produced by the Ban 
D2 line, which has been clearly inter preted in term s of the HFS shown by two of the 
seven stable isotopes of barium (see IStenflol 1 19971: iBelluzzi et al. ][1qo3). Given that 
the only stable isotope of scandium shows HFS, and noticing the similarities between 
the two lines in the intensity spectrum, it was believ ed that th i s sign al could also be 
interpreted in terms of HFSq However, as shown by IBelluzzil (2009), this is not the 
case: indeed the theory unexpectedly fails to model this signal, even though the effects 
of HFS are correctly taken into account, and the same modeling assumptions used for 
barium are considered. Such a failure has importa nt consequences : if the Sc 11 signal 
is produced by a physical mechanism neglected by IBelluzzil ((2009), or that the theory 
does not account for (for example PRD effects), then not only it will be important to 
identify this mechanism, but it will be also interesting to understand why it seems not 
to be necessary in the case of barium. Perhaps the problem on the physical origin of the 
triplet-peak structure of the Ban signal is still waiting for a complete answer! 



4. Theoretical approaches accounting for PRD effects 

The most natural way to describe PRD effects is by means of a redistribution function, 
a powerful tool which naturally appears in th e "traditional " scattering approach to the 
physics of resonance polarization (see, e.g., Stenflo 1994). This approach, which has 
been successfully applied for treating PRD phenomena, has, on the other hand, the 
limitation of being able to describe only two-level atoms with unpolarized, infinitely 
sharp lower level. The redistribution function H is defined so that 

-> -> . dfl\ dQi 

K(vi,Q2,* 2 ;vi,Qi,*i)dvidvi— --r - (8) 



3 According to this classification the "M" signals are all those signals which show positive polarizing peaks 
(or lobes) in the wings of the line, a depolarization approaching to the line center, and (eventually) a sharp 
polarizing peak in the line core. 

4 It should be observed that only 18% of the barium isotopes and 100% of scandium show HFS. This is 
in perfect agreement with the fact that the central peak of the barium signal is much higher than the other 
peaks, while the three peaks of the scandium signal have the same amplitude. 
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represents the probability that in a scattering process an incoming photon of frequency 
v\, propagation direction Q.\, and polarization state e\ is scattered in a photon with 
frequency v%, propagation direction Q.2, and polarization state £i- If the scattering of 
a radiation beam, whose polarization properties are described through the 4 Stokes 
parameters, is considered, then the redistribution function is replaced by a 4 x 4 redis- 

tribution matrix 7?j,(v2, "2! Vi, ^i). 

When the frequency dependence can be ignored, the redistribution matrix reduces 
to the so-called scat t ering phase m atrix (or Rayleigh phase matrix) Pij{0.z', (see, 
e.g.. iHamiltonll 19471 : IStenfloll 19941 . or LL04). In the atomic rest frame the angular and 
frequency dependencies can be separated, so that the redistribution matrix can be factor- 
ized into the product of the scattering phase matrix times a scalar function describing 
the frequency redistribution. This is no more true in the laboratory frame where the 
Doppler effect, due to the motion of the atoms, has to be taken into account. In this 
case, the scalar function describing the frequency redistribution also depends on the 
directions of the incoming and scattered bea ms, and an intri cate coupling between fre- 
quency and angle dependencies appears (see lHummerlll962l for the unpolarized case). 
In order t o simplify the proble m, the so-called hybrid approximation has often been 
used (e.g. jRees & Saliball 1982h . This approximation, which is justified only by heuris- 
tic arguments, consists in calculating the average over the angles of the scalar function 
describing the frequency redistribution, so that the angle and frequency dependencies 
can be factorized as in the atomic rest frame. 

From the 1970s, the effects of PRD in frequency have been investigated using 
different approximate forms of the redistribution matr ix, genera l ly bas ed on the angle 
averaged type I and typ e II redistribution functions of iHumirierl ( 1962 ) (see Sect. 1 of 
ISampoorna et al.l |2010L for a his torical review). Particularly useful results were ob- 
tained by iFaurobertl (11987L Il988b who investigated the role of PRD effects considering 
both angle-averaged and angle-dependent redistribution functions, also including the 
effect of a magnetic field. These investigations pointed out the limitations of approx- 
imate forms of the redistribution matrix previously proposed, and showed that the ap- 
proximation of complete redistribution (CRD) is generally adequate to describe the line 
core polarization. The fact that PRD effects (i.e. coherent scattering) are important in 
the wings of strong resonance lines, while they are usually neglig ible in the l ine core 
has been then confirmed by several other investigations (see, e.g.. lFrischlll996b . 

From the theoretical point of view, the most tangled physical aspect in a scatter- 
ing process of polarized radiation is the role of collisions (elasti c and inelast i c). A 
first, fundamental investigation on this aspect was carried out by lOmont et al.1 ( 1972 ) 
who int roduced the role of collisions in th e quantum theory of Ram an scattering devel- 
oped b y lFiutak & Van Kranendonkl dl962l) . Starting from this work. iDomke & Hubeny 
1988)' derived a very general analytic expression of the redistribution matrix for scat- 



tering of a rbitrarily polarized radia tion by an atom undergoing collisions. Using this 
expression iFaurobert-SchollI (1 1992b point ed out the sensitiv ity to the elastic collisional 
rate of the line-wing polarization, while iNagendral (119941) showed how depolarizing 
collisions affect the line-core polarization. A detailed investigation on the role of var- 
ious atmospheric parameters on the emergent line ar polarization, c a lculate d includ- 
ing PRD effects, has been recently carri ed out by Sam poorna et al.l (l2010l) . To this 
end, ISampoorna & Trujillo Buenol (1201 dl) developed new, very efficient numerical ap- 
proaches for PRD radiative transfer based on the Gauss-Seidel and Successive-Overrelaxation 
iterative methods. 
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Since the frequency redistribution matrix does not explicitly appear in the statis- 
tical equilibrium equations and in the radiative transfer equations, in order to take into 
account PRD effects within the theoretical approach described in Sect. 3, it is necessary 
to consider the next terms in the perturbative expansion of the atom-radiation inter- 
action (i.e. second and higher order processes). Unfortunately, such a generalization 
of the theory has revealed to involve many difficulties, since the equations that come 
out when higher order terms are included become immediately extremely complex and 
practically unmanageable, unless further approximations are introduced. 

Nevertheless, first important results have already been obtained. Under the sim- 
plifying assumptio n of a two-level a tom with unpolarized lower level, and neglecting 
stimulation effects, B ommierl (1 1 997a!) was able to perform all the calculations including 



higher order terms in Eq. (|7), thus deriving a generalized expression for the radiative 
transfer coefficients. With respect to the expression obtained when only first-order pro- 
cesses are taken into account (valid under the complete redistribution approximation), 
the new emission coefficient contains an extra term (the so-called Rayleigh scattering 
term) which descri bes the contribution of co herent scattering. Remarkably, the redis- 
tribution matrix of Domke & Hubeny (1988) is recovered if the generalized emission 
coefficient derived by iBommierl (11997 ah is expressed in terms of a redistribution matrix: 

e,<V2,4)= j/°>-§ I dvi (h— ±Yj R ij(y2>ti2;vi,&i)Sj(y u tii) 

V 4 

+7 /0) -feB P (vo)0(vo-v) £,•,(), (9) 
v o 

with i - 0, 1, 2, 3 (corresponding respectively to the Stokes parameters /, Q, U, and V), 
?7 (0) the line integrated absorption coefficient, vo the line frequency, and Sj (J = 0, • • • , 3) 
the four Stokes parameters. The last term describes the contribution to the emitted 
radiation coming from collisionally excited atoms, the branching ratio being e - r ^+r R ' 
with Tj and F R the inelastic collisions de-excitation rate and the radiative de-excitation 
rate, respectively. The 5/ t o shows that this contribution is unpolarized (collisions are 
assumed to be isotropic). The quantity Bp represents the Planck function in the Wien 
limit, while <ft (vo - v) is a norm alized Lorentzian profile. The redistribution matrix is 
given by (see [b ommierj l 1 997bh 

2 

Rij = 2 W K (J f , J u ) [P {K \&2, do] [a r n (v 2 , vi) + [/3 {K) - a] r m (y 2 , vi)) . (10) 

This expression is valid in the atomic rest frame, where, as can be observed, the fre- 
quency and angular dependencies can be factori zed. The quantity Wk is a numerical 
factor introduced by Landi DeglTnnocentil (Il984l) . the matrix P^ K ^ is t he A~-th multipole 
compo nent of the Rayleigh phase scattering matrix, also introduced by lLandi DeglTnnocentil 



( 1984). The functions rjr(v2, vi) and r///(v2, vi) are given by 

r//(v 2 ,vi) = 0(v o - vi) 8(y 2 - Vi) , (11) 
rm(v 2 ,v{) = 0(v o - vi) (f>(v 2 - vi) , (12) 



where (/>(vq - vi) is a Lorentzian absorption profile, while 5{v2 - v\) is the Dirac delta. 
The first function describes the coherent contribution to the scattering, the second one 



12 



Luca Belluzzi 



the contribution in complete redistribution. The branching ratios are given by a and 
\0K) _ a ^ respectively, with 

a= T . and ff*> = -ft — , (13) 

where T^ is the elastic collisions rate (responsible for line-broadening and for destruc- 
tion of frequency correlations between incoming and scattered photons), and D (K ^ is the 
A'-multipole depolarizing rate due to elastic collisions (responsible for depolarization). 
The quantity a is the probability that a radiative decay takes place before any kind of 
collision: it thus represents the branching ratio for the scattering processes that are co- 
herent in the atom rest frame. The quantity 0- K \ on the other hand, is the probability 
that a radiative decay occurs before an elastic depolarizing collision. The difference 
_ a ] thus gives the probability that the radiative decay takes place after the atom 
undergoes an elastic collision that redistributes the photon frequency, but that does not 
destroy atomic polarization: it thus represents the branching ratio for the scattering pro- 
cesses that are completely redistributed in frequency in the atom rest frame. It should 
be observed that a branching ratio term (1 — e), taking into account the fraction of 
atoms that are radiatively excited, is already included in the branching ratios a and 
[jg(*0 _ a j (f or a disc ussion on the physical m eaning of the various quantities appearing 
in Eq. ( fTfJl) , see also Sampoorna et al.ll2Q10h . In this context, the CRD limit is obtained 
when the elastic collisions are so efficient that any frequency correlation between the 
incoming and outgoing photons is completely destroyed. Indeed, when Te»Tr the 
coherent term disappears (a — > 0), and the resulting redistribution matrix corresponds 
to the one which can be obtained within the theoretical framework described in Sect. 3 
(valid under the flat spectrum approximation). As already observed in Sect. 3, when the 
elastic collisional rate is veiy high, depolarizing collisions are also very efficient, and 
polarization phenomena become negligible. For this reason, when the CRD limit due 
to high collisional rates is considered, it is customary to assume the rates Fe and D^ K) 
as independent, though they have the same physical origin. On the other hand, in the 
limit of no collisions (Te = = 0) the branching ratio \ft K ^ - a] vanishes, and the 
scattering is completely cohere nt, as implied by e nergy conservation. 



Within the same approach, iBommier] (Il997bf) also investigated the effect of a mag- 



netic field of arbitrary strength (which is found to mix angle and frequency dependen- 
cies also in the atom rest frame), as well as the transition to the laboratory frame. 

Another interesting theoretical approach, based on the density matrix formalism, 
which can be applied to multi-level atoms, in the presence of arbitrary magnetic fields, 
and that doe s not require the flat spectrum a pproximation to be assumed, has been 



proposed by ILandi DeglTnnocenti et al.l (119971) . In this approach the density-matrix 



formalism is generalized to the framework of the so-called metalevels theory, in which 
the atomic levels are assumed to be composed by a continuous distribution of infi nitely 
sharp sublevels. The formulation presented in lLandi DeglTnnocenti et all (1 19971) does 
not take collisions into account, a nd thus describes the limit of pure coherent scattering. 
This scheme has been applied by lLandi DeglTnn ocentil d 19981) for the interpretation of 
the peculiar signal observed in the second solar spectrum of the sodium doublet. 

More recently, "intermediate" theoretical schemes, exploiting the strength points 
both of the scattering approach, and of the density matrix formalism have been pro- 
posed. An interesting example is the approach proposed by iHolzreuter et al.l (120051) for 
the investigation of the physical origin of the triplet-peak structure shown by several 
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signals of the second solar spectrum. Although these approaches are able to account for 
many fundamental physical aspects, in general they are not derived in a self-consistent 
way from basic physical principles. This makes extremely important, and in some 
cases also very difficult, to clearly establish their limits of applicability, which might 
often result to be fairly restrictive. 

5. Conclusions 

While until the end of the 1990s, the various branches of solar polarimetry, from the 
instrumentation to the observational techniques, from the theory of polarization to its 
solar application, were developed hand in hand, in the last ten years, as it can also 
be noticed from the references mentioned in this paper, progresses in the field of pure 
theory are proceeding significantly slower than in the other research branches. 

In 1995, when the first "Solar Polarization Workshop" was organized, the need of 
developing theoretical tools suitable for describing the effects of frequency redistribu- 
tion was already clear, and the problem of deriving a general analytical form of the re- 
distribution function was strongly debated. Today, fifteen years later, our understanding 
of this physical aspect has substantially improved, different promising approaches have 
been proposed, and for the case of a two-level atom (with unpolarized lower level) the 
PRD problem seems to be solved. Nevertheless, the development of a self-consistent, 
general theory for radiative transfer in the presence of PRD effects (applicable to multi- 
levels atoms, and accounting for the fundamental role of quantum coherences) remains 
since many years a fundamental unsolved problem. We are aware of the difficulties that 
the formulation of such a theory implies, and it is reasonable to expect that its devel- 
opment will proceed rather slowly, through a series of small, demanding steps. The 
strong possibility of not achieving novel and outstanding results in short or medium 
periods might explain why not many people (in particular among the new generations) 
decide to cope with this difficult challenge, definitely necessary for a new substantial 
acceleration in this research field. 

The theory developed by Landi DeglTnnocenti at the beginning of the 1980s, and 
explained in great detail in LL04, keeps representing the most solid, advanced and 
complete theoretical approach to the physics of spectral line polarization developed so 
far. More than twenty years after its first development, after a countless number of 
successful applications, its soundness for the interpretation of all those signals which 
can be treated within the limit of complete frequency redistribution appears to be well 
established. 
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